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Abstract 

We use the "dispersive method" to investigate non-perturbative effects in the 
transverse momentum distribution of vector bosons produced in pp collisions. 
The assumption of a non-perturbative modification to the running coupling at 
low scales leads to additional contributions in impact parameter space propor- 
tional to —b 2 logQ 2 and — b 2 . Our results, which we expect to be valid provided 
t is not close to 1, are shown to account for a substantial proportion of the 
total non-perturbative contribution extracted from data. 
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1 Introduction 



The measured transverse momentum distributions of W and Z bosons at the Teva- 
tron and LHC will provide great insights into the structure of hadrons. Indeed some very 
interesting data has already been released 0. On the theoretical side, the application 
of perturbation theory has led to a full next-to-leading order calculation and next-to- 
leading log resummations [§, |], [5], |j . These are merged to give what is at present our best 
perturbative prediction for the distribution. There are two different approaches in the lit- 
erature fU to the resummation of these logarithms: either the contributions are resummed 
directly in transverse momentum space, or a Fourier transform is performed and the calcu- 
lation developed in impact parameter space. These approaches were compared in [|j] and 
shown to differ only for NNN-leading logs. In this article we will use the impact parameter 
formalism, since it is more readily suited to the calculation at hand. Perturbation theory 
has however proved insufficient to describe experimental data, and it is necessary to include 
in calculations large non-perturbative effects. Parametrisations of these contributions have 
been proposed ||, and values for parameters have been extracted from data 0, |§; however 
a full theoretical understanding of these effects has yet to emerge. 

One technique that has had success in describing non-perturbative effects over a wide 
range of QCD observables is the dispersive approach M. It is related to the renormalon 
model (see |10| for a review) but rests on slightly different foundations. The primary 
assumption is that one can define a universal running strong coupling a s (k 2 ) that has no 
singularities in the complex plane except for a branch cut along the negative real axis. In 
particular there is no Landau pole, which is an unphysical artefact of perturbation theory. 
The contribution already accounted for in a fixed-order perturbative calculation is then 
subtracted off, leaving what is termed the 'non-perturbative' correction. The assumption 
of universality of the running coupling, and thus of the non-perturbative parameters derived 
from it, is most easily tested using data from event shapes in e + e~ annihilation [[ll|] and 
deep inelastic scattering |JT2|: these are particularly large corrections all proportional to the 
first moment of the non-perturbative modification to the coupling. Studies indicate that 
universality is approximately realised, both giving support to the model, and prompting 



potential refinements of some calculations [ 13 



Power corrections to observables in the Drell-Yan process have previously been stud- 
ied using the renormalon and dispersive approaches in |14], [15|, [16[ . In addition, non- 
perturbative corrections to a calculation involving a next-to-leading log resummation were 
studied for the first time in the dispersive approach for the e + e~ energy-energy correlation 



T7fl . Here we apply the same technique to the next-to-leading log resummation of the W/Z 
transverse momentum distribution. A simplified overview of the perturbative resummation 
is presented in section [|, while the power corrections to the distribution are calculated in 
section 131. There then follows in section HI a discussion of the results. 



2 Transverse momentum distribution 

We recall the perturbative prediction for the transverse momentum distribution in 
two stages: first we consider the emission of a single gluon; and then we resum leading 
logarithms arising from multiple collinear gluon radiation, figure [I]. Non-perturbative cor- 
rections from the modification to the coupling at low scales can then be added to this 
perturbative result. 
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Figure 1: Vector boson production with collinear gluon emission. 



2.1 Single gluon emission 

Suppose in the collision of two hadrons P and P', a quark and antiquark annihilate into 
a vector boson. The parton model gives the total cross-section for the process as 

ff ( r ) = H I dX<iX ' [<l( X > /^V, tf) + V( X , V 2 )} Vqq(x, X\ T, fi 2 ) , (2.1) 



where q and q are the annihilating species of quark and antiquark respectively, whose 
distribution functions in P and P' appear above, a qg is the parton-level cross-section, and 
the standard Drell-Yan variable is 

( f , (2-2) 



T = 



where Q is the boson mass and ^/s the total energy in centre-of-momentum frame. We 
have also explicitly introduced the factorisation scale fi 2 , which is in principle arbitrary but 
which (it will be seen below) is naturally chosen for the transverse momentum distribution 
to be (b /b) 2 . 

In order to find the differential p_i_ distribution da/d 2 p±, where p_i_ is the two-component 
transverse momentum of the produced boson, normal to the incident hadron momenta, we 
simply insert a factor S 2 (p± + J2i k_u) into the phase-space integral, where here % runs over 
all emitted gluons. It is convenient to move into impact parameter space, writing 



«S 2 (px + i:k lj )=/(05e«><P i+ E, ki ,). 



(2.3) 



Thus we find 

da , N r d 2 b 



d 2 p, 



■)-/ 



3 ib-p ± 



V / dxdx \q(x, [i 2 )q'(x, /i 2 ) + (<->)] a q g(x, x', r, /i 2 , b) , (2.4) 
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where a qq is now the cross-section calculated with an additional factor e lb k± for each real 
emitted gluon. Performing the angular integrals then gives 



da 



(2.5) 



roo i>l 

/ dbbJ (bp±)^2 dxdx' q(x, p?)q {x' , p 2 ) + (<->) a q g(x, x', r, p 2 , b) . 



At Born level, O(a s ), we find 



a (0) = 

V an 



qq 2xx's N r 



2ty5(xx's - Q 2 ) ATraqgQ' 



47T 2 « 



qq 



N r s 



5(xx' — t) , 



(2.6) 



where 4ira qq is the coupling of the quark q and antiquark q to the vector boson, which in 
the case of W production includes the relevant CKM matrix elements. 

The 0(a s ) correction comprises contributions from real and virtual gluon emission. 
Denoting the 4-momenta of the incoming quark and antiquark p and p', that of the gluon 
k and the vector boson q, we may write in the centre-of-momentum frame of the incoming 
quark and antiqark, 



P = -^-r— (1,0,0,1) 



p' = ^-r— 1,0,0,-1 



(2.7) 



and 



k— (k°, k± cos 0, k± sin 0, A; 3 ) . (2.8) 
To evaluate the real contribution, we note that the squared matrix element is given by 



jbkj_ I 



M | 2 = e ib -^-- [Ana qq ] {Ma s C F ) 



N, 



(p ■ q) 2 + (p' ■ gf 
(p ■ k) {jpf ■ k) 



(2.9) 



which, in terms of the variables 



z = 



XX 



k 2 

Q 2 ' 



(2.10) 



may be written 



ib.k X | M ,2 = ibQVAco S <j> 4?r 8-K 2 a qg a s C F 
11 zA N r 2tt 



2 + 2z 2 - 4zA 



(2.11) 



The Lorentz-invariant integration measure and phase space for the squared matrix elements 
is 

9(1 - z ) /-(i-*)V(4*) zdA 

(z.lzj 



47T JO 

and thus we obtain the 0(a s ) real contribution as 
„ rcal _ 4^-0(1 - z) a s C F f(^)V^) dA 



l-z) 2 - AzA J o 2vr 



qq 



N r s xx 1 



2tt 



/ 

Jo 



Jfhn /r> U^-M 
A \/(l-z) 2 -AzA 



(2.13) 
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This is divergent at small A when the gluon becomes soft or collinear with one of the 
incoming quarks; we will therefore find it convenient to separate this logarithmic divergence 
from the finite remainder. 

In order to simplify the calculation we take the Mellin transform 

E(N,pl) = ( X dTT^n^vl) ■ (2-14) 
Jo 

This yields 

1 roo AlT^O - 

nN.vl) = - dbbJoibp^^r^bNQ'N + QNQ^^ + Cii^b) - R(N,b)} , (2.15) 

/ JO qq JVcS 

where in the final term on the right hand side, the 1 represents the Born level contribution, 
and the quantity C\ — R the finite and divergent parts of the 0(a s ) correction. The 
moments of the parton distributions are 

q N = f 1 dxx N ~ 1 q{x) . (2.16) 



Thus we obtain for the real gluon emission: 



n Ct? f°° dA . /•1+2A-2 A /A(1+A) 9 i 9 7 2 _ a A 

Cr l _ R re,l = /• dJ\j {bQ ^ f V ^ zN _, 2 + 2z AzIX 

2tt Jo A 'Jo . I{\ _ z y _ AzA 



(2.17) 

The z integral may be evaluated as a series expansion in A (see the appendix for details) 
to give the collinear and soft divergences: 

^ r ° al (^6) = ^/ o 1 ^Jo(&gv^){21ogA + 45 1 (iV) + ^ + ^} . (2.18) 

Returning now to the virtual contribution, by expressing the loop integrals as a single 
integral over we obtain 



*« ~ N c s 6{XX T) 2vr Jo A Ke \ 3+ v^4A l0g v^4A + 1 / " ^ 

This also diverges at the lower limit, and thus for the purposes of a full 0(a s ) calculation 
a regulator is required, such as dimensional regularisation or, as performed in section ^|, a 
small gluon mass. But this expression is sufficient for obtaining the logarithmic divergences. 
The virtual contribution is separated into finite and divergent pieces according to: 

47T 2 (T - r l 

aZf = ^-f^Sixx' - t) C7 irt - R vlIt ] , (2.20) 
N c s L J 

where the collinear divergence R is found by expanding the integrand in A: 

wirt a sC F fidA 



R v 



jf — {-21ogA-3}. (2.21) 



271 

Thus the sum of real and virtual divergences becomes: 



(2.22) 
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There is a certain freedom in the choice of R, notably in the choice of the upper limit 
kj_ = Q 2 . But this is compensated by the finite contribution C±, which is defined such that 
C\ — R gives the full 0(a s ) virtual contribution. 

The radiator R contains the logarithmic divergences associated with single-gluon emis- 
sion. It is a non-trivial function of N, but can be simplified by a wise choice of factorisation 
scale. We begin by splitting the radiator into three pieces 



where 



R = R 1 + R 2 + R 3 , (2.23) 



) 2 AV2 



a s C F fQ dk\ ^ .„ 2 

2tT Jo k'j L ~ i - N N 



R2 = -=^r L -TT^-M^iS-AS^N)---^^) (2.25) 



a s C F fQ 2 dk 2 , ( „ , , 2 2 \ . 

The first two pieces may be simplified by making the de facto approximation of the 
Bessel function: 

[1 - J {bk±)] « B(bk ± - 6 ) , (2.27) 
where 60 = 2e -TB and 7^ is Euler's constant, as described more fully in the appendix. Thus 

I* ^ (3 _ 4S[(iv) _ 2 2 — \ 

2tt y ( 6 /fe)2 fci V N + l) K ' 

The contribution i?3 is formally divergent. Using dimensional regularisation, with 4 + 2r] 
dimensions, it becomes 

a s C F [ Q2 dklU-K^/k 2 )" 1 ! / 2 2 \ . 

R ^^rl # r(i + ;i ( 4 ^) + ^ + FTT- 3 + H < 2 - 30 > 

where /1 is the factorisation scale which appears in the quark distributions. On application 
of Ms, we obtain 

where the finite contribution arising from the 0{rj) term in equation ( |2.3(J| ) has been ab- 
sorbed into C\. Let us therefore choose the factorisation scale to be /i 2 = (b /b) 2 , in order 
that the terms R2 and R3 cancel. Alternatively, since these contributions are proportional 
to the non-singlet anomalous dimension, starting from any factorisation scale /i 2 the terms 
i?2 and -R3 may be absorbed into the parton densities by means of a scale change to (bo/b) 2 . 
Thus the single-gluon contribution may be written as 

m>p±) = 0/ dbbMbp ± )j:-^[ qN ((bo/bf)U(bo/b) 2 )+(^ 

x {l + Ct{N,b)-R{b)} , (2.32) 
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where the collinear and soft logarithms are represented by the radiator 

and the remaining 0(a s ) contribution is to be found in C±. 

2.2 Multiple gluon emission 

Consider now the case in which many gluons are emitted before the vector boson is 
produced. Each splitting contributes to the total transverse momentum of the event: we 
are interested only in the contributions from soft and collinear logarithms, which may be 
calculated using the collinear approximation. 

Let the 4-momenta of the incident hadrons be P and P' . In the centre of momentum 
frame we may write 

P = ^(1, 0, 0, 1) P' = ^(1, 0, 0, -1) , (2.34) 

where s is the total energy of the collision. 

Suppose the (anti)quark in P has momentum fraction x; and let it radiate n gluons 
with momenta ki, (i — 1, • • • , n). We may write the Sudakov decomposition 

h = aiP + piP' + k ±i , (2.35) 

in terms of which the integration measure is 

(2 - 36) 

Let Zi be the fraction of remaining longitudinal momentum retained by the (anti)quark 
at the ith splitting, i.e. 

on = (1-Zi) (x-J2aA (2.37) 



x — J2j=i a j 



Zi = - J:}' . (2.38) 



Then the collinear-divergent pieces of the squared matrix element are given by the usual 
splitting function: 



\M\U = ^r^P^^-\M\h^ m *, (2-39) 
2tt J 2tt Zi 1 — Zi 



where 0j is the azimuthal angle of the ith gluon. Since we are considering collinear emission, 
the phase space is subject to the condition > fa, which translates in our approximation 
to 

*<1-|. (2.40) 
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This boundary is not the boundary of physical phase space; rather it is imposed as part 
of the collinear approximation. It is necessary therefore to check that any dependence 
in our calculation on this limit is a genuine physical effect and not an artefact of the 
approximation. 

We may treat radiation from the other incoming (anti) quark in the same way, but with 
aj and /?< interchanged. The real part of the partonic cross-section appeaing in Q2.5|) is 
then: 



real 1 1 [A 1^ 1 TT f OL s Cp f dk 2 ±i f dz { 1 + z\ ■ 



n'=0 - i'=l 

Taking the Mellin transform leads to: 



/■oo 4-7T rv - 

S(iV,pi) = - / d66J (6p±) E + QnQn] e~ Rm) , (2.42) 

I Jo iV r S 



99 ' 

where the real emission radiator is 



R Ies i(N, b) = / -prMok^) dzz - . (2.43) 



7T JO kj_ JO 1 — Z 

We then take into account virtual gluon radiation using the universal virtual contribu- 



tion to the radiator 17 



/IN a s C F [Q 2 dk\ /-i-*±/<3 , l + z 2 , 
RvirtuM = / -pr / dz- , (2.44) 



1 JO 

so that the full radiator is 

R(N, b) = H^r f a ' d 4 r tl ' Q i'— V - «"-'JM»x)l • (2.45) 

7T JO ft^ JO 1 — Z 1 J 

Performing the z integral at small k\ yields: 

R(N, b) = ^ f ^G(k 2 JQ 2 , bk ± , N) . (2.46) 

where 

G(k 2 JQ 2 ,bk ± ,N) = -2[l-J (bk ± )]lo g %-3 + 2J (bk ± )(2S 1 (N) + ±- + ' 



Q2 V iV JV+1 

Q W 



+4 [1 - NJ (bk ± )) ^- + (%) . (2.47) 



It should be noted that: 



This result is strictly valid provided N is not too large, which here means approx- 
imately N < Q/k±, since as N — > oo the contribution from real emission in (|2.45| ) 
vanishes. But we are not interested in such large values of N, since they correspond 
to values of r close to 1, while for vector boson production at the LHC or Tevatron 
the order of magnitude is r ~ 0.01. 
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The terms of order k\jQ 2 do not generate logarithmic divergences: they are thus 
included in the finite contribution to the cross-section and therefore do not feature in 
the radiator. The term of order k±/Q is a spurious artefact of the collinear approxi- 
mation, that is not present in the single-gluon result. The collinear approximation is 
thus valid only for logarithm resummation, and may not be used for determinations 



of power corrections |15 



The decomposition of the radiator proceeds just as in the single-gluon case. The fac- 
torisation scale is naturally chosen to be /i 2 = (bo/b) 2 and we obtain the exponentiated 
form: 



1 r°° 

2 Jo 



dbbMbp^^r 11 

1 V r 5 



QN((bo/bY)q' N ((bo/bY) + 



(2.48) 



where the leading collinear and soft logarithms are represented by the radiator 



R 



a,Cj, 



dk 2 



-2*4-3 



(2.49) 



2tt J{bo/Vf k\ 

and the remaining contribution to fixed order in perturbation theory is given by the factor 

C{a s ) = 1 + d(N,b) + 0{a 2 s ) . (2.50) 

Next-to- leading logarithms have been calculated || || ||, § , and may be included in the 
radiator by writing 



dk 2 



2vr J(b /b) 2 kj 



, k 2 N 
- 21 ° g #- 3 , 



(2.51) 



where the perturbative coupling is defined in the Bremsstrahlung scheme and is, to 
next-to-leading order, 



4tt 



where 



aT(Q 2 ) + [aT(Q 2 

(67 - 3it 2 )C A - 207^71/ 



Q 2 
k 2 



— lo S TT + TZ 



FT 



K 



9 



(2.52) 



(2.53) 



Having now recalled the perturbative contribution to the transverse momentum distri- 
bution, we are in a position to calculate non-perturbative effects from the modification to 
the coupling at low scales. 



3 Non-perturbative contribution 

The dispersive approach to power corrections has been applied to a variety of QCD 
observables. The basic postulate is that we may define a universal running strong coupling 
a s (k 2 ) that is well-behaved in the infra-red as well as in the ultra-violet, having no singu- 
larities in the complex plane apart from a branch cut along the negative real axis. We may 
then write a dispersion relation 

(3-D 
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where p s {p 2 ) is the discontinuity of a s across its branch cut: 

p s (/i 2 ) = ^K(eV)-«s(e"V)} 



(3.2) 



Non-perturbative contributions to the transverse momentum distribution arise from the 
modification to the running coupling at low scales: 



a s (k 2 ) = a^(k 2 )+5a s (k 2 ) 



(3.3) 



where the contribution af T , given in ( |2.52| ), is that already accounted for in the perturbative 
part of the calculation. In the dispersive approach the corresponding modification to the 
spectral function, 5p s (k 2 ), gives the correction to the coefficient function 



SC{N, b) 



2tt 



00 dk 2 
~k 



2 5p s (k 2 ) C(k 2 /Q 2 )-C(0) 



(3.4) 



where (a s CF/2Tr)C(k 2 /Q 2 ) is just the single-gluon contribution to C(N,b) calculated as 
though gluon has mass k 2 . As customary we define e = k 2 /Q 2 . 

The coefficient function has contributions from both real and virtual gluons. Consider 
first the real contribution. A calculation of the massive single-gluon matrix element and 
phase space gives 



roo rl/\ , fZ(A,e) 

c-i = r ^ Jo(6Q va) dz 

Jo A Jo 



_! 2A(A + z 2 A - 2zA 2 + 2ez 2 + 2e 2 z 2 + e 2 z 2 A) 



(A + ez) 2 J(l -z + ez) 2 - 4z(e + A) 



where the upper limit on the z integral is 

(l + 2A + e )-2^/( e + A)(l + A) 

z(A ' e) = o^? 



(3.5) 



(3.6) 



Note how this reduces to perturbative case (|2.17| ) when we put e = 0. 

In order to proceed with this integral we expand the Bessel function. We require only 
those terms that are non-analytic in e as e — > 0, and these are generated at the lower limit 
of the A integral. So the expansion is legitimate for our purposes and we obtain 



C 



real / 



l/(l+7i) 2 



dz z 



N-l 







{l-z+ezf/(Az)~e ( b 2 Q 2 A 

— + 



X 



2(A + z 2 A - 2zA 2 + 2ez 2 + 2e 2 z 2 + e 2 z 2 A) 



(A + ez) 2 J(l -z+ ez) 2 - Az(e + A) 



(3.7) 



In order to evaluate these integrals we follow the method of 0. The integral over A is 
relatively straightforward, giving 



C 



real / 



1/(1+^) 2 



dz z 



N-l 



1 — z — ez 



tanh 1 



1 — z — ez) 2 — Aez 2 



z — ez 



Q(z,e)J(l -z-ez) 2 - Aez 2 



(3i 
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where the functions P and Q are given by 

P(z,e) = A + A(l + e) 2 z 2 + 2eb 2 Q 2 z(l- z + z 2 -ez + 2ez 2 + e 2 z 2 ) (3.9) 
b 2 2 

Q(z,e) = 4+ -^-(l + z + z 2 + ez + 10ez 2 + e 2 z 2 ) . (3.10) 
bz 

We then make the substitution 



v (1 - z - ez) 2 - 4e^ 2 1 

r] = tanh -2 z = — — — — , (3.11) 

1 — z — ez 1 + e + lyje cosh 77 

which leads to the expression for the characteristic function 

C real (e) = AI N + 2tb 2 Q 2 I N+l + 4(1 + e) 2 I N+2 - 2e(l + t)b 2 Q 2 I N+2 + 2e(l + efb 2 Q 2 I N+3 
-\tb 2 Q 2 J N+1 - 16eJ N+2 - §e(l + e)b 2 Q 2 J N+2 - §e(l + lOe + e> 2 Q 2 , (3.12) 

in terms of the integrals 

77 tanh 77 



Ijv(e) = / drj 



M<) = L d V Tl , : , - jJ—^ . (3.13) 



1 + e + 2 x Acoshr/) 7V 

sinh 2 rj 
(1 + e + 2 x /ecoshr]) J 

Using the values of these integrals given in the appendix, we obtain the contribution from 
real gluon emission as: 

C real (e) = {vr 2 + (25 1 (Ar)+loge) 2 -4S 2 (iV) + 2(l + ] ^- T )(25 1 (iV)+loge)} 
b 2 Q 2 ( 1 1 



Ntt z 



6 ViV + 2 N-l 



+e {2(2S' 1 (JV) + loge) (n 2 + 3iV - 3 + 

-N(2Si{N) + log e) 2 + 4NS 2 (N) - 6N 2 -6N + 16- ^qjy } 



C,rQ {vr 2 + (2S' 1 (JV) + log e) 2 - 45 2 (iV) + 2 - U 



4 [ ■ v ^ / ■ o / -v j N + 2 

+2 (| + - l) (2ft(JV) + loge)} + • • • . (3.14) 

Now we turn to the virtual contribution. Since the virtual gluon does not contribute 
to the transverse momentum of the vector boson, we obtain the same as that calculated in 
for the total Drell-Yan cross-section, which is 

C virt = (1 + e) 2| 2Li2( _ e)+21ogelog(1 + e) + ^_ log 2 e |_7_ 2e _ (3 + 2e)loge 

= (-log 2 e-31oge + y-^ + e (^-2 log 2 e + ^ - 4^) + 0(e 2 ) . (3.15) 

The power correction is most conveniently expressed in terms of moments of the modi- 
fication to the coupling, 5a s (k 2 ): 

SC(N,b) = i^l -^6a s (k 2 )V(k 2 /Q 2 ) , (3.16) 
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where T>(e) is the discontinuity of C(e) across its branch cut: 

1 



2ni 



C(ee l7T ) - C(ee 



(3.17) 



The leading term in the characteristic function is the loge divergence proportional to the 
anomalous dimension: this is subtracted off and absorbed into the parton density functions. 
From the finite contribution we obtain 



Vie) 



-2e | (iV + 2) log e + 2iV5i (N) — N 2 — 3N + 3 
eb 2 Q 2 



1 



loge + 25' 1 (JV) + -|+ 1 



N N + 2 



N + l 



(3.18) 



The b- independent term on the first line gives the 1/Q 2 corrections to the total cross-section 
calculated in ||. Of particular relevance to the transverse momentum distribution is the 
6-dependent term on second line, which gives the correction 



5C(N,b) (logQ 2 - 2S X (N) - ^ - 



1 - A 



N + 2 

The non-perturbative parameters A 2 and A' 2 are defined in terms of 5a s by: 



(3.19) 



Z7l J fh I 



The extent of our knowledge of A 2 and «A' 2 is (I7|, [LJ| that *4. 2 
present unconstrained but appears consistent with zero. 



(3.20) 
(3.21) 

0.2 GeV 2 and A is at 



4 Results and discussion 



In the previous section we derived the non-perturbative 1-gluon contribution to the 
transverse momentum distribution in impact parameter space, which by performing the 
inverse Mellin transform may be written: 



1 



2 Jo 



dbbJ {b P± )e- R{b) K{T,b) [l - b 2 W(r,b) 
where the perturbative contribution is 

dxdx' q(x, (bo/b) 2 )^ (x', (b /b) 2 ) + (<->) 8(xx' — r) 



K(T,b)=J2 4n2aq " ' 



N s 

qq x ^c° 



and the non-perturbative component is 



W(t, b) = i {A 2 log Q 2 + A 2 - A 2 } - ^-K(r/z, b) (l 



+ z z 



(4.1) 



(4.2) 



l-z) 



(4.3) 
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Figure [| shows a plot of W(t, b) against r at fixed values of b, for on-shell Z° production 
with A<i = 0.2 GeV 2 , A' 2 = 0, evaluated using the MRST (central gluon) distributions [pOf] ; 
the slight dependence of W(t, b) on b is logarithmic, coming only from the scale of the 
parton distributions. It is also clear from the figure that the result given here is not valid 
near r = 1: it becomes negative and diverges logarithmically. This is because our expression 



for the Mellin transform was valid for N not too large (equation (|2.47f) in the perturbative 
case, and analogously though not so obviously in the non-perturbative calculation), i.e. 
N < Q/k±, while as r — > 1 the large moments become important. The approximation is 
valid away from this region, and is certainly legitimate for LHC or Tevatron data. 




Figure 2: Graph of W(r, b) against r for b = 1 GeV 1 (solid line), b = 0.1 GeV 1 (dotted) 
and b = 0.01 GeV -1 (dashed), evaluated at the Z mass with A 2 = 0.2 GeV 2 and A' 2 = 0. 



We thus predict a non-perturbative correction that behaves quadratically in b near 
b = 0, in agreement with existing models based on a Gaussian suppression factor ||. We 
can make our result Gaussian also by exponentiating: 



1 -b 2 W(r,b) « e 



-b 2 W(r,b) 



(4.4) 



Whether or not the power correction does exponentiate is an open question — intuitively 
it seems reasonable, since we expect contributions from all soft gluons, not just the gluon 
nearest the hard vertex; yet it is not proven and should be treated with care, since exponen- 
tiation normally applies to collinear and soft logarithms, while the power corrections come 
from terms that converge in the infra-red. Nevertheless there are several good reasons for 



at least studying the exponentiated correction. Firstly, the integral (|4.1|) diverges at large b 
due to the non-perturbative contribution growing as b 2 , while the exponentiated expression 
is convergent. At small b, where the non-perturbative correction is small, whether or not 
we choose to exponentiate makes minimal difference. 

Secondly, we can consider the difference between the exponentiated and un-exponenti- 
ated expressions to give a measure of the theoretical uncertainty in our calculation, since 
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we expect higher orders to contribute at approximately this level. Thus at small b the error 
is very small, while at large b it becomes huge. This has the same pattern as the general 
impact parameter formalism, which suffers from considerable uncertainties at large b due 
to the parton distributions becoming unknown. 

Finally, there have been some fairly reasonable phenomenological fits of Gaussian mod- 
els to data, such as performed in ||. A two-parameter fit yielded the non-perturbative 
contribution to the exponent 



-b 2 



0.34 GeV 2 log 



Q 



0.16 GeV^ 



.1 GeV, 

while a three-parameter fit, with an additional linear term, may be written 

Q 



0.48 GeV 2 log 



0.41 GeV 2 



+ 6 [0.09 GeVlog(100a;a; / )] 



,1 GeV, 

These may be compared with the prediction from the dispersive approach, which is 

Q 



b z W 



0.2 GeV 2 log 



1 GeV 



0.2 GeV' 



at 



0.01 



(4.5) 



(4.6) 



(4.7) 



where we have used A 2 = 0.2 GeV 2 and A' 2 = 0. It is seen that this accounts for a 
substantial portion of the total non-perturbative function fitted to the data; however the 
r-dependence of the fitted forms ( [4.5|) and (|4.6|) differs so greatly from that predicted here 
that it is hard to give more than this qualitative statement. A full fit to data is beyond the 
scope of this current article, but such a fit would be a further test of the dispersive model of 
power corrections, would enable one to extract a value for the unknown non-perturbative 
parameter A' 2 , and would test the assumption of universality which predicts A2 = 0.2 GeV 2 



from DIS structure functions 18 



An alternative approach to the non-perturbative contribution ( |3.19| ) would be to inter- 
pret part of it as a shift in the scale at which the parton distributions are evaluated. We 
may write the contribution as 



5C(N, b) 



b 2 A 
2C F 



iV + 2 



where the non-singlet anomalous dimension is 

'l + z 2 \ 



dz z 



N-l 



C f 



2S l (N) - 



1 



1 



(4i 



N N + 1 



(4.9) 



The first term on the right hand side of ( |4.8| ) then shifts the factorisation scale from (bo/b) 2 
to (b /b) 2 + 5fi 2 , where 



b 2 A, 
2C F 



(bo/b) 2 



A 1 



2tt 



(4.10) 



If as a first approximation we replace the integrand by its value at the lower limit, we 
obtain the shift 



b 2 A 2 \ a s ((b /b) 2 )C f 
2tt 



-3 GeV 2 



(4.11) 
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where in obtaining the approximate numerical value, the running of a s is neglected. 

So a negative correction to the factorisation scale is predicted, in addition to which 
there is the remainder of the non-perturbative contribution calculated above. Thus instead 
of ( FOp , we may write 



1 r°° -> r 

£(r,pi) = 77 / dbbJ (bp ± )e- R ^K(r,b) 1 - b 2 W(r,b) 

2 Jo 1 



(4.12) 



where the perturbative contribution is now 

~q(x, (b /b) 2 + 5/2 2 )q'(x', (bo/b) 2 + 5fi 2 ) + 6(xx' - r) , 



HI 



N r s Jo 



and the non-perturbative component is 



(4.13) 



W{t, b) = - {2A 2 log Q 2 -A 2 - 2A' 2 \ - C -K(r/z, b)(l-z + 2z 2 ) . (4.14) 

4 > 2K(r, b) Jt z v ' 

Figure ^| shows a plot of W(t, b) against r at fixed values of b. Again there is a slight 
logarithmic dependence of W(t, b) on b, coming from the scale of the parton distributions. 
Note however that W as a function of r is very much flatter than W was, since we have 
absorbed the most singular behaviour into the parton density functions. However this 
advantage is to some extent offset by the fact that the shift of about -3 GeV 2 in the 
factorisation scale reduces further the region of impact parameter space that may be treated 
perturbatively. 



W 




J I I L 



0.4 0.5 
X 



Figure 3: Graph of W(t, b) against r for b = 0.5 GeV" 1 (solid line), b = 0.1 GeV" 1 (dotted) 
and b = 0.01 GeV -1 (dashed), evaluated at the Z mass with A 2 = 0.2 GeV 2 and A' 2 = 0. 
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Appendix 



In this appendix are included some technical details of the calculations, 
(i) Integral (gT?| ): 



We require a series expansion in A for the integral 
L(A) = 



l+2A-2y / A(l+A) 



dz z 



v-1 2 + 2z 2 - AzA 



1 - z) 2 -AzA 



(A.l) 



To this end we divide the integration region at z — 1 — y/A/n, where vA C k < 1, 
and take each part separately. Since k is an arbitrary quantity, we consistently discard 
all terms that vanish as k — >■ 0. 



For < 2 < 1 - vA/k we may expand the integrand directly in A to obtain 



Lx(A) 



1-VA/k 



dz z 



N-l 



2z 2 



.aiog^-^AO-l- 



+ 0(A) 
2 



K 



N N + l +W — + 0(A )' (A ' 2 » 



where 



N-l 



n=l 



rv 



(A.3) 



The contribution from the second integration region can be evaluated by substituting 



z = 1 + 2A — 2-y/A(l + u + A) and expanding in A to obtain 

2 



L 2 (A) 



^+v/A) -l-A du 
(I Ju 



Vl + u 



ANVA + 0(A) 



A 



-2 log KT - AN— + 0(A) . 



(A.4) 



Adding together the two contributions ( |A.2|) and (|A.4j) gives the infra-red divergent 
part as shown in ( |2. 18| ) . Note also that the term in a/A cancels — the fact that 
in the collinear approximation a 1/Q term remains indicates that it is not a valid 
approximation for studying power-suppressed corrections. 
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(ii) Approximation ( |2.27|) : 

The contributions Ri and i?2 to the radiator are of the generic form 

R *-Jo ~k\ 
Let us define the function 

rQ 2 r/A-7 
'/i 2 



[1 - Jo(bk ± )} fi(k] 



F^ 2 ,Q 2 ) = - [[ 

J a- 1 



(A.5) 



(A.6) 



Then, integrating ( A.5|) by parts and making the substitution x = bkj_, we see that 

rbQ 



R; 



dxFi((x/by,Q 2 )J' (x) 



(A.7) 



Further, by inserting the appropriate functional form of f%(kj_) into (|A.6| ), we see that 
Fi((x/b) 2 , Q 2 ) is some function of log(x/b) which is slowly varying in x compared with 
the oscillating Jq(x). We may therefore extend the range of integration to infinity, 
and expand in logx about x = b : 



F^x/b) 2 , Q 2 ) = F^bv/b) 2 , Q 2 ) + F{((b /b) 2 , Q 2 ) \og(x/b ) + ■ ■ 



(A. 



where here F[ = dFi/ dlogx, and the higher derivatives have correspondingly fewer 
logarithms. Then, since 



rod 

/ dx J'q(x) = —1 

J 

we obtain from ( |A.6| ) 



Ri 



Q 2 dk 2 



dx Jq(x) log(x/& ) = 



(A.9) 



(A.10) 



/(6 /6) 2 k\ 

Thus the Bessel function 1 — Jo{bk±) may be replaced in the radiator by the step 
function Q(bk± — bo) to the accuracy we are using. 



(iii) Integrals ( |3.13|) : 

By substituting u = cosh 77, we may write 

r°° du log(w + \/u 2 — 1) 1 



1 u (l + e + 2^u) N 

A 

du- 



dz \og 2 (z+ Vz^l) 



Ani Jc z (1 + e — 1\flz\ 



N 



Me) h --(l + e + 2Veu)^ 
where the contour C is: 



2m Jc 



dz 



1 + e - 2jlz\ 



N 



,(a.h; 




>Re 
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This contour encloses poles in the integrand at z = and z = (1 + e)/(2y / e). The 
residues of these poles may be evaluated by expanding around them, yielding the 
results: 



1 



1 - iVe) 7T 2 - 4S 2 (iV) + (2S 1 (N) + loge) 



3iV 2 — N — 2 — (N — 1)(N + 2)(2S 1 (N) + loge) + 0(e 



1 

Te 



(N — 1)(N — 2) 



+ e(25 , i(JV)+loge-l 



N 



+ 0(e 2 ) 



.(A.12) 
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